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Abstract 

Power corrections both to the strong coupling constant and to the structure 
function itself are estimated on the basis of the LO, NLO and NNLO QCD 
analysis of xF% structure function data . The sign of correction to the coupling 
constant is found to be negative. The x-shape of a higher twist contribution to 
the structure function is stable up to NNLO. 

PACS numbers: 12.38.Bx; 12.38.-t; 13.85.Hd ;14.20.Dh 



Power corrections to perturbation QCD predictions both for Q 2 -evolution of the 
running coupling constant and the structure function itself have intensively discussed 
recently 0, 0. At present, precise measurements of structure functions (SF) and de- 
tailed theoretical calculations of QCD predictions for scaling violations provide an 
important means of accurate comparison of QCD with experiment. On this basis, 
during the last years the important information on the x-shape of power corrections 
to the structure function of a nucleon was obtained by the QCD analysis of data on 
deep- inelastic scattering of leptons || - || . 

In the present note, detailed LO, NLO and NNLO QCD analyses of precise xF 3 
structure function data have been carried out in order to evaluate the power cor- 
rection to the QCD running coupling constant 

*s{Q 2 ) = *f cD m+%// (i) 
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where in NNLO the coupling constant a p s (Q 2 ) can be expressed in terms of inverse 
powers of L = \n(Q 2 /Aj^) as 

vQCD (n 2\ 1 (3iln(L) 

as {Q) = M"wr (2) 

+ 0^\ftln 2 (L) - d\ln{L) + (3 2 (3 - f3 2 } 

Notice that /3b = 11 - 0.6667/, ft = 102 - 12.6667/, (3 2 = 1428.50 - 279.611/ + 
6.01852/ 2 . 

We perform the QCD analysis using the method of the Jacobi polynomial expansion 
of structure functions. This method of solution of the DGLAP equation was proposed 
in H and developed both for unpolarized |9|] and polarized cases [[10] . The main formula 
of this method allows approximate reconstruction of the structure function through a 
finite number of Mellin moments of the structure function 

h(r) Nmax n , s , s 

xF»<~ (x,Q 2 ) = ^ + x a (l- X y £ e a n >P(x)J2cf(a,/3)M j+2 (Q 2 ) (3) 

^ n=0 j=0 

The Q 2 -evolution of M N {Q 2 ) is defined by the perturbative QCD 



d N 

H N (QlQ 2 )M% CD (Ql), N = 2, 3,... (4) 

d N = 7 (0) ' 7V /2/3o 

Here a s (Q 2 ) is the strong interaction constant, >y^' N are nonsinglet leading order 
anomalous dimensions, and the factor Hn (Qq,Q 2 ) contains next and next-to-next 
to leading order QCD corrections. Power corrections to the coupling constant are 
introduced formally in accordance with (JJ). 

Unknown coefficients D {Qq) in (ED could be parametrized as the Mellin mo- 
ments of some function: 

M® CD (N,Q 2 ) = f 1 dxx N - 2 ax\l-x) c , iV = 2,3,... (5) 
Jo 

The shape of the function h(x) as well as parameters A 2 , a, b, c, and are found 
by fitting the experimental data on the xF 3 (x, Q 2 ) structure function J7|. Detailed 
description of the fitting procedure could be found in [|J. Both terms h(x)/Q 2 and 
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A2/Q 2 are considered as pure phenomeno logical. For a possible analytic expression see 
|L1| . The target mass corrections are taken into account to the order o(M^ ud /Q 4 ). 
The results of the fit are presented in Table 1 and Figure 1. 

The values of A 2 in Table 1 are opposite in sign to the lattice results for NLO and 
NNLO @] A% LO = 0.22(2) GeV 2 and A^ NLO = 0.21(2) GeV 2 obtained in different 
renormalization schemes. Notices also that in |12| "lattice data" for the coupling con- 
stant were analysed, whereas in the QCD analysis of structure functions the anomalous 
dimensions and coefficient functions were involved too. Even in the leading order, the 
first factor in the right-hand side of (||) reproduces the terms with powers differing 
from —2, which initially appears in ([!]). However, the absolute value of the parameter 
A 2 and a large value of the scale parameter A.j^g are in qualitative agreement with []T2j . 



On the other hand, the negative value of A 2 is in agreement with predictions of [111 





A 2 [GeV 2 ] 


Am^ [MeV] 


x 2 




h(x) 


- free 




LO 


-0.261 ± 0.053 


1140 ± 110 


70.5 / 96 


NLO 


-0.130 ± 0.027 


788 ± 92 


71.7 / 96 


NNLO 


-0.123 ± 0.017 


561 ± 46 


73.4 / 96 




h(x 


)=o 




LO 


-0.137 ± 0.011 


834 ± 31 


115.6 / 96 


NLO 


-0.049 ± 0.012 


584 ± 69 


116.5 / 96 


NNLO 


-0.046 ± 0.013 


561 ± 74 


103.2 / 96 


NLO 


-0.011 ± 0.008 


267 ± 36 


135.4 / 96 


NNLO 


-0.023 ± 0.005 


290 ± 36 


125.5 / 96 



Table 1. The results of the LO, NLO (N Max = 10) 
and NNLO (N Max = 6) QCD fit (with TMC) of 
xF 3 data 0. (Qg = 3 GeV 2 , Q 2 > 3 GeV 2 , f=4). 
The bottom two lines correspond to substitution 
of (U) into the moments of the coefficient function 
only. 

The values of constants A 2 for NLO and NNLO are approximately the same in 
agreement with the statement that the 1/Q 2 corrections to all orders in a s are of the 
same order [] [0 . 

1 The LO result should not be considered on this matter because a nontrivial contribution of the 
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The shape of h(x) slightly differs from the results of analysis in Ref. || with A 2 = 0. 
The effect of decreasing the power correction to the structure function |5], |6| while going 
from LO to NNLO of the perturbative QCD does not exist, as canbeseen from Fig.l. 

A special fit for the case h(x) = gives a negative value for the parameter A 2 . 
The increase of the x 2 parameter shows that using only power corrections to coupling 
constant, one could not reach a good description of experimental data and it is nec- 
essary to introduce additional power corrections to the structure function itself. A 
large difference between the value of A2 for h(x) = and h(x) 7^ indicates strong 
correlations between two power terms: A2/Q 2 and h(x)/Q 2 . 

Even higher x 2 is obtained when the formal substitution of ([!]) is applied to co- 
efficient finctions and is not applied to the anomalous-dimention-dependent factor of 
structure function moment. This result is presented in the bottom two lines of Table 1. 
The values of A 2 are small and negative, but parameter the A^- is in agreement with 

In conclusion it should be noted that the presented values of A 2 should be taken with 
a great caution. One should make use od strict analitic expressions for the Q 2 evolution 
of the structure function moments with power corrections to coupling constant for a 
reliable fit of data. The nuclear and threshold effects should be taken into account, as 
well. 
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Fig.l The results of the LO, NLO and NNLO extractions of twist-4 contributions 
of h(x). The power correction to the QCD running coupling constant as(Q 2 ) = 
af CD {Q 2 ) + A 2 /Q 2 is included. 
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